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We consider the thermodynamic consistency of the charge response function in the (extended) 

Hubbard model. In DMFT, thermodynamic consistency is preserved. We prove that the static, 
homogeneous DMFT susceptibility is consistent as long as vertex corrections obtained from the 
two-particle impurity correlation function are included. In presence of a nonlocal interaction, the 
problem may be treated within extended DMFT (EDMFT), or its diagrammatic extension, the dual 
boson approach. We show that here, maintaining thermodynamic consistency requires knowledge of 
three- and four-particle impurity correlation functions, which are typically neglected. Nevertheless, 
the dual boson approximation to the response is remarkably close to consistency. This holds even 
when two-particle vertex corrections are neglected. EDMFT is consistent only in the strongly 
correlated regime and near half-hlling, where the physics is predominantly local. 

PACS numbers: 71.10.-w,71.10.Fd,71.30.-|-h 


I. INTRODUCTION 

Dynamical mean-field theory (DMFT)^^^ has emerged 
as an important tool for the description of strongly corre¬ 
lated electron systems. The theory treats the local tem¬ 
poral correlations and can be applied at arbitrary interac¬ 
tion strength, allowing it to capture the Mott transition. 
Extended DMFT (EDMFT)^^® was developed to treat 
nonlocal interactions, that can lead to phenomena such 
as charge-ordering. Due to recent technical advances, 
EDMFT has received renewed attention®^^^ and has been 
adapted to spin systems (so-called spin-DMFT)^^. At 
the same time, diagrammatic extensions of DMFT and 
EDMFT have appeared, such as the dynamical vertex ap¬ 
proximation (DTA)^^, dual fermion^^d^^ one-particle 
irreducible (IPI) approach^^, EDMFT-|-GW^dO,i8^ dual 
boson (DB)^® and TRILEX^*’. These are advanced, nev¬ 
ertheless approximate methods to the correlated many- 
body problem. For any approximation it is desirable 
to fulfill certain requirements, such as the conserva¬ 
tion laws of energy and charge^^d2 thermodynamic 
consistency^^d2^ While these cannot always be fulfilled 
in practice, it is nevertheless important to know for the 
interpretation of the results, if and how well these re¬ 
quirements are fulfilled. In this work, we address the 
thermodynamic consistency of the linear charge response 
in DMFT, EDMFT and DB. 

The response to an external perturbation is relevant 
for many experiments and for understanding the under¬ 
lying physics. The linear response is related to corre¬ 
lation functions of the system via the celebrated Kubo 
formula^^. The correlation function is often more at¬ 
tractive to theorists. Especially for space- and time- 
dependent external perturbations (e.g. fluctuating mag¬ 
netic fields), the direct response may be difficult to cal¬ 
culate, because it involves the non-homogenous and non¬ 
equilibrium properties, whereas correlation functions of 


the unperturbed homogeneous equilibrium system are 
more accessible. In this case the Kubo formula provides 
a connection between the two and between theory and 
experiment. The correspondence between the response 
and correlation functions is a property of the exact solu¬ 
tion of the system and this equality is not automatically 
satisfied by approximate solutions^^^^®. In fact, this kind 
of consideration was very important in the development 
of the theory of the electron gas^^. 

We call the linear response of an approximation ther¬ 
modynamically consistent if the two ways, the direct cal¬ 
culation by varying the external field and that from the 
correlation function, yield the same result. The so-called 
zero- and one-particle quantities of DMFT are known 
to be thermodynamically consistent^®’^®, i.e., the same 
density is obtained from deriving the grand potential 
and from the Green’s function. The DMFT magnetic re- 
sponse®° of the double-exchange model also turns out to 
be thermodynamically consistent®^, while EDMFT and 
the related spin-DMFT®® appear to violate thermody¬ 
namic consistency. When both the response and the cor¬ 
relation function can be calculated, thermodynamic con¬ 
sistency can be checked. 

Here we apply such analysis to the charge susceptibility 
of an (extended) Hubbard model with local and nonlo¬ 
cal density-density interaction. We find that the charge 
response in DMFT can be obtained in a consistent man¬ 
ner provided vertex corrections from the impurity two- 
particle correlation function are included. In EDMFT 
and DB, thermodynamic consistency is in general vio¬ 
lated and we show that the knowledge of the three- and 
four-particle impurity correlation functions would be re¬ 
quired to restore this deficiency. While this is often im¬ 
practical, the DB approximation turns out to contain the 
dominant contributions to yield a response that is consis¬ 
tent to very good approximation, over a broad interaction 
and filling range. We identify those contributions. 
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The paper is organized as follows: The model and the 
Kubo formula for the charge response are introduced in 
Sec. 11. In Sec. Ill, we examine the case of local inter¬ 
action and show, by resorting to the DB approach, that 
DMFT is thermodynamically consistent in the Hubbard 
model. We then use thermodynamic consistency as a 
yardstick to measure how other approximations to the 
charge correlation function, namely DB without vertex 
corrections, approximation by the impurity susceptibil¬ 
ity and the DMFT bubble, perform. In Sec. IV, we ap¬ 
ply the same approach to the extended Hubbard model 
and find that three- and four-particle interactions beyond 
a two-particle ladder summation are necessary to reach 
thermodynamic consistency. Finally we summarize our 
findings in Sec. V. Additional technical facets and a de¬ 
tailed derivation of the DMFT response formula can be 
found in the Appendix. 


II. MODEL 

We study a single-band, two-dimensional extended 
Hubbard model on a square lattice away from half-filling. 
The Hamiltonian of the extended Hubbard model is given 

by 

Uk)(T j 

+ '^Unj^nji + ^'^Vjknjnk, (I) 

3 (jk) 

with Uj = rija = ^ '1' 

k denote site indices and c|-^ (^o-) creation (an¬ 

nihilation) operators for an electron with spin a on site 
j. In this work, we restrict ourselves to nearest-neighbor 
hopping and interaction, i.e., tjk = t and Vjk = V when 
j and k are nearest-neighbors and zero otherwise. The 
physical parameters are the hopping integral t, local in¬ 
teraction U, nonlocal interaction V and chemical poten¬ 
tial /i. For our calculations, we take the half-bandwidth 
D = At as the energy unit and use inverse temperature 
P = 20. As in EDMFT-|-GW^®, it is convenient for DB 
calculations to write the nonlocal interaction in Eq. (1) 
in terms of density fluctuations, which requires a Hartree 
shift M ™ ^be chemical potential. 

We study the dependence of the density (n) on the 
chemical potential using the linear response formalism. 
The Kubo formula relates the density response to a 
two-particle correlation function. The electron com¬ 
pressibility d{n) ld\i is related to the static homoge¬ 
neous correlation function as (Appendix A): d(n) /d/i = 
- limq^o lim^^^o 


III. LOCAL INTERACTION 

We first address the thermodynamic consistency for 
the simpler case of the Hubbard model (V = 0) treated 
within DMFT^. In DMFT, the lattice problem is mapped 
to an auxiliary single-site impurity problem with action 

*^imp [c , c] — ^ ^ /r 

I^tT 

( 2 ) 

CD 

Here v and a; are the fermionic and bosonic Matsub- 
ara frequencies respectively, = 'Yliv and a 

normalization factor /3 is implied in the sums over fre¬ 
quencies. The hybridization function describes the 
electron hopping processes from and to the impurity and 
includes the effect of the lattice in a mean-field manner. 
It is determined self-consistently via the condition 

^ ^ 9 1^(7-J (1^) 

k 

where is the Green’s function of the impurity and 
denotes averaging over momenta. The lattice Green’s 
function is given in momentum space as 

^kJcr = 9ua + ^I'o- - tk, (4) 

where tk is the Fourier transform of the hopping. 

A. Density response 

In DMFT, the impurity density is the same as the lat¬ 
tice density as a direct consequence of the self-consistency 
condition. The impurity and lattice compressibilities, 
however, are different, since the impurity density de¬ 
pends on the chemical potential /r directly, but also via 
the hybridization function which is determined self- 
consistently. As a consequence, the compressibility de¬ 
composes into the impurity compressibility and a part 
originating from the variation of A^o- due to a change in 
chemical potential: 

d{n) _ a(n)™P a(n)™P 

dfj. ~ dn ^ ^ dA^^ 

Here and in the following we use (.) and (.)™^ to denote 
lattice and impurity averages, respectively. The right- 
hand side consists of three parts. The first is the impurity 
compressibility. The second part is proportional to a two- 
particle impurity correlation function d{nY'^'^ /dA^a oc 
- + (n)™P (c*^c^^)™P. According to 

Eqs. (3) and (4), the information about the lattice comes 
in through the variation of the self-consistent hybridiza¬ 
tion, i.e., dA^^jd^i. By working out the details, as done 
in Appendix B, the DMFT response can be expressed 
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Figure 1: (Color online) Density (n) obtained in DMFT as a 
function of chemical potential /r for different values of the local 
interaction 17, at t = 0.25, T = 1/20. The Hartree contribu¬ 
tion jj,H is subtracted from the chemical so that half-hlling 
occurs at /r — fin = 0 for all interaction strengths. Close to 
the Mott transition {U = 2.35) the chemical potential be¬ 
comes flat near half-filling due to the development of a gap. 


in terms of one- and two-particle impurity correlation 
functions. A crucial aspect in the derivation is that the 
impurity problem is solved (numerically) exactly using 
continuous-time quantum Monte Carlo^^“^^. The im¬ 
purity response and correlation functions are therefore 
thermodynamically consistent. The final result is the 
static, homogeneous susceptibility of the DB approxima¬ 
tion, which, for the Hubbard model, has been proven^^ 
to be exactly equivalent to the DMFT susceptibility in¬ 
cluding vertex corrections^®. 

Before confirming this numerically, we show the de¬ 
pendence of the electron density (n) on the chemical po¬ 
tential /r at various interaction strengths U in Fig. 1. 
At U = 0, the exact density can be obtained directly 
from the non-interacting density of states. At finite U 
the DMFT densities are shown. The curves approach 
each other in the dilute limit. For low density, the sys¬ 
tem’s properties, including (n), depend weakly on U be¬ 
cause the physics is determined mainly through the ki¬ 
netic term. The chemical potential is shown with the 
Hartree contribution = U{n) /2 subtracted, so that 
half-filling occurs at ^ — fin — 0 for all values of U, due 
to the particle-hole symmetry of the bipartite square lat¬ 
tice. The DMFT correlation function is evaluated at a 
fixed value of p, to find the compressibility (orange circles 
in Fig. 2). On the other hand, the DMFT compressibil¬ 
ity can also be obtained by numerically computing the 
derivative of the density with respect to fj, of the curves 
in Fig. 1. Here we calculate it as the difference quotient 
A{n)/Afj,, which introduces a negligible error, as can be 
judged from the plots. It is depicted by the blue squares 
in Fig. 2. 

Let us focus first on these data. As a general observa¬ 
tion, Fig. 2 shows that the electron compressibility of the 
weakly and strongly interacting systems behaves differ- 
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(c) U = 2.35, V = 0 

Figure 2: (Color online) Compressibility d{n}/d/j. for the 
Hubbard model {V = 0), computed within various approx¬ 
imations and three different values of the on-site interaction. 
The chemical potential is shown at the bottom of each graph, 
the corresponding density (see Fig. 1) is shown at the top. 
Blue squares denote the response function computed by vary¬ 
ing the chemical potential as An/Ap, while other results are 
computed from the correlation function in different approx¬ 
imations. Note that only (n) < 1 (less than half-filled) is 
shown, results for (n) > 1 can be obtained by particle-hole 
symmetry. 
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ently as half-filling is approached. The results are con¬ 
sistent with Fermi liquid theory, i.e. the compressibil¬ 
ity decreases with increasing interaction and decreasing 
spectral weight at the Fermi level. At U = 2.35 (cf. 
Fig. 2c), the compressibility is correspondingly low. In 
Fig. 1, the U = 2.35 graph is indeed almost flat near half¬ 
filling. Note that the compressibility does not diverge as 
half-filling is approached^®’^^, because the temperature 
is above the critical temperature of the metal-insulator 
transition. 

Comparing the results for the difference quotient with 
DMFT/DB, we see that also numerically, the thermo¬ 
dynamic consistency is apparent at all values of U as 
expected (we recall that the DB result is the same as 
DMFT for = 0). 

We further show results for the compressibility ob¬ 
tained from the correlation function approximated by 
a bubble of renormalized DMFT propagators (blue di¬ 
amonds in Fig. 2). For small interaction, the approxi¬ 
mation performs well, as expected. Already at moder¬ 
ate values of the interaction however, the result deviates 
strongly. In particular, it does not describe the small 
compressibility close to half-filling. Irrespective of the 
value of the interaction, the result approaches the one of 
DMFT in the dilute limit, where vertex corrections are 
expected to be small. It is further instructive to deter¬ 
mine the compressibility from the local impurity suscep¬ 
tibility. This approximation neglects the fact that the im¬ 
purity density depends on the hybridization, cf. Eq. (5). 
While it contains local vertex corrections, it neglects non¬ 
local ones and the momentum dependence contained in 
the bubble. As a result it performs poorly at weak inter¬ 
action. For moderate and strong interaction, it is close 
to the response function only near half-filling, where the 
physics is predominantly local. 

In Fig. 2 we examine another approximation, which 
in a sense interpolates between the results of the bub¬ 
ble and impurity susceptibility. This approximation has 
been introduced in Ref. 35. Without going into all de¬ 
tails here, we mention that it is obtained by neglecting 
the two-particle impurity vertex and the corresponding 
ladder diagrams contributing to the susceptibility in the 
DB expression for the susceptibility [see also Appendix B, 
in particular Eq. (B13)]. It however still includes a three- 
leg vertex and the bubble diagram, which is renormalized 
through this vertex. In terms of DB perturbation theory, 
this diagram contributing to the dual bosonic self-energy 
is second order in the three-leg vertices (dual electron- 
boson interaction) and hence has been coined “second- 
order DB”. We provide some more details on this ap¬ 
proximation in Sec. IV. 

The results of the calculation are shown as red tri¬ 
angles in Fig. 2. At small U, the fermion-fermion in¬ 
teraction is small and neglecting the two-particle vertex 
works well. Furthermore, at U = 0.5 the three-leg ver¬ 
tex is essentially independent of fermionic frequency^®, 
so that second-order DB is expected to behave similarly 
as the bubble. It agrees almost perfectly with the differ¬ 


ence quotient. At large interaction strength, the agree¬ 
ment remains remarkably close. Close to half-filling, we 
have seen that the local susceptibility captures the es¬ 
sential physics. However for quarter filling and below, 
we see that both vertex corrections and the momentum 
dependence of the bubble are important. Neither the 
local susceptibility nor the bubble alone reproduce the 
charge response. From the derivation of the DMFT sus¬ 
ceptibility in Appendix B, we see that the three-leg ver¬ 
tex emerges from the variation of the impurity density 
with respect to the hybridization d (n) jdA^„ [Eq. (B5)] 
and the variation of the impurity Green’s function with 
respect to the chemical potential dgi,„/d^ [Eq. (B 6 )]. 
These terms are apparently important beyond the bub¬ 
ble in the correlated regime. On the other hand, ne¬ 
glecting the two-particle vertex corresponds to neglecting 
dgvaldAyiai [Eq. (B7)], i.e. the change of the impurity 
Green’s function with respect to The agreement is 

hence only approximate. Such an approximation has a 
further shortcoming, as it violates the Ward identity and 
charge conservation®®. We note however, that this ap¬ 
proximation is computationally signihcantly less expen¬ 
sive than the full DMFT correlation function, because the 
two-particle impurity vertex function 7 does not have to 
be computed. In this particular case, the three-leg vertex 
is determined as a function of fermionic frequency only, 
at a single bosonic frequency uJn=o- This is comparable 
in computational cost to DMFT calculations. 

Finally we note that the above discussion is relevant for 
the case of a local self-energy. Consistency will in general 
be violated if the self-energy is momentum dependent (as 
in dual fermion), see Ref. 39 for an example. To maintain 
consistency in that case requires additional diagrams in 
the formula for the correlation function. 


IV. NONLOCAL INTERACTION 

In this section, we examine thermodynamic consis¬ 
tency in the extended Hubbard model within EDMFT 
and DB. In EDMFT, the effect of screening due to the 
nonlocal interaction is accounted for through a retarded 
interaction in the impurity model. This frequency depen¬ 
dent interaction is determined through a self-consistency 
condition analogous to the one for Ay^r in DMFT. We 
recall that the resulting impurity action reads 

5'imp[c*,c] = + g. - Ay„]cy^ 

U(7 

+ U'^ ^ X! (6) 

UJ LO 

and the self-consistency conditions are 
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(a) U = 0.5, V = 0.2 
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Figure 4: Diagrammatic representation of nq'f?o!tj='o • 

the self-consistency conditions (7), (8) are fulfilled. 

The DB approximation can be viewed as a diagram¬ 
matic extension of EDMFT. It allows us to incorporate 
additional diagrams (polarization corrections^'^) into the 
charge susceptibility. To this end, (9) is replaced by 

^quj ~ [XuJ T Xui^quiXuj] T (1^) 


(n) 


C\J<D->-h>-C\J00-;J-O(DC\J05L0CMa3'<^-0^'^C0 

oo->-->-c\Jc\JcO'^'^LOLn(Dr^h>.coooo>05 

dddddddddddddddddd 



(b) U = 1.5, V = 0.2 


( 72 ) 


C'J<D-i-(D->-h>.<NC0-^0><DC'JC0'^0>'^05<NLr)C00>O 

OO-i-->-C\JC\JC0C0-^'^L0«D(Dh>.r^C0C00^05050^O 

dddddddddddddddddddddi-^ 



(c) U = 2.35, V = 0.2 


Figure 3: (Color online) Compressibility as a function of 
chemical potential and density for finite retarded interaction. 
Labels are as in Fig. 2. 


Here Xbj denotes the impurity susceptibility and 

Xq^=XZ^+^u7-Vq (9) 

is the EDMFT charge correlation function. In our calcu¬ 
lations, we solve the impurity model (6) repeatedly until 


where Ilqi^ is dual bosonic self-energy. Note that 
EDMFT is recovered for flq^ = 0. We evaluate B di- 
agrammatically. The leading diagram to Ilq;^ is shown 
in Fig. 4. It is second order in the electron-boson interac¬ 
tion and we refer to this approximation as second-order 
DB approximation. The ladder DB is obtained by replac¬ 
ing one of the triangles through a renormalized triangular 
vertex containing ladder diagrams^'*. The explicit expres¬ 
sions are given by Eq. (B14) and Eq. (C2) respectively. 

Although diagrammatic corrections to the self-energy 
and more elaborate self-consistency schemes are possible 
in DB^o, we restrict ourselves to “one-shot”-type calcu¬ 
lations here (see also the discussion in Sec. V). I.e., we 
solve the EDMFT equations, and only after convergence 
evaluate the susceptibility according to (10) instead of 
(9). As a result, the impurity quantities and in particu¬ 
lar the Green’s function and density (and hence An/A/i) 
are the same in EDMFT and DB. 

In EDMFT and DB, we can determine the variation of 
the density with respect to the chemical potential analo¬ 
gously to the foregoing. We have 


d{n) 

dfi 


dfi ^ dfi 

^ dA^ dfi ■ 


( 11 ) 


An additional term contributing to the compressibility 
emerges due to the retarded interaction A^^. As dis¬ 
cussed in the previous section, the contributions in the 
first line are accounted for in DB, but the one due to the 
retarded interaction is not. Performing steps similar to 
the derivation in Appendix B shows that cl(n)™'’ /i9A^ oc 
— (n)™*’ ■ In fact, there will 

even be a four-particle correlator in the derivation, aris¬ 
ing from dxui/dA^i in the variation of the self-consistency 
condition (see Appendix D). Hence, in spite of the fact 
that the Kubo formula relates the response to a two- 
particle (lattice) correlation function, three- and four- 
particle impurity correlation functions appear in the ex¬ 
pression for the response. This is due to the presence of 
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the retarded interaction. For most applications, account¬ 
ing for three-particle or even higher-order correlations is 
impractical. It is therefore important to examine their 
contribution to observables. There has been some dis¬ 
cussion on the impact of three-particle interactions in the 
literature^^’^^, but in numerical studies they are usually 
neglected^^’^®’^® (with the exception of Ref. 42 where the 
effect was found to be small). 

Because of the higher-order correlation functions in the 
response, EDMFT and DB cannot fulfill thermodynamic 
consistency exactly. In order to examine to what extent 
they are consistent, we show the electron compressibility 
for the extended Hubbard model for different local and 
fixed nonlocal interaction strengths in Fig. 3. We see 
that the overall compressibility is smaller for finite V at 
a given U (cf. Fig. 2). The agreement with A (n) /A/r 
is remarkably close for DB for all values of U. The ef¬ 
fect of three- and four-particle correlations appears to 
be negligible, even close to the Mott transition. The 
strongest deviations occur for relatively low filling, where 
they remain small (cf. insets). As before, the EDMFT 
susceptibility is consistent in the correlated regime close 
to half-filling, where the physics is predominantly local. 
In this regime, we have found corrections of DB pertur¬ 
bation theory to EDMFT to be small^°. We have also 
performed a calculation at significantly lower tempera¬ 
ture of /? = 50. Without nearest-neighbor interaction, 
this interaction strength U = 2.35 would be very close 
to the metal-insulator transition. We do not find any 
qualitative difference to the foregoing. 


V. CONCLUSION AND DISCUSSION 

We have assessed the thermodynamic consistency of 
the charge response in the (extended) Hubbard model 
within DMFT, EDMFT and DB. We have proven for the 
Hubbard model, that DMFT yields a consistent response. 
The DMFT bubble is consistent only at small interac¬ 
tion. On the other hand, the local impurity susceptibil¬ 
ity is consistent when the physics is predominantly local, 
i.e., at large on-site interaction and close to half-filling. 
Beyond these two regimes, additional diagrams are re¬ 
quired to yield a consistent charge response. We found 
that dominant contributions are retained even when ne¬ 
glecting nonlocal vertex corrections from the two-particle 
vertex (this amounts to a second-order approximation in 
terms of dual variables). From these results we conclude 
that the most important ingredients for maintaining a 
consistent response is the momentum dependence at low 
interaction, local vertex corrections at strong interaction 
and their interplay at intermediate interaction, whereby 
a renormalized electron-boson coupling given by the im¬ 
purity three-leg vertex is essential. 

In the extended Hubbard model with nearest-neighbor 
interaction, we have shown that the response is deter¬ 
mined not only by two- but also by three- and four- 
particle correlations. This is a consequence of the self- 


consistent retarded interaction included in the impurity 
model. The DB approximation, which neglects associ¬ 
ated diagrams, however, yields a response which is con¬ 
sistent to good approximation. Also here, non-local long- 
range vertex corrections appear to have a small effect. At 
least in this aspect, the effect of third- and fourth-order 
correlation functions is negligible, which is important in 
practice. EDMFT is, in general, not consistent. It is con¬ 
sistent when the physics is essentially local. We found 
that the most important contribution beyond DMFT to 
correct this deficit is a second-order diagram in the renor¬ 
malized triangular vertex, which effectively includes mo¬ 
mentum dependence into the EDMFT polarization. 

The above observations are valid in the case when no 
diagrammatic corrections to the fermionic self-energy are 
taken into account, i.e. it remains local and equal to the 
impurity self-energy. For a non-local self-energy, main¬ 
taining consistency requires additional contributions to 
the correlation function. 
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Appendix A: Linear response formalism 

The particle number and the density response can be 
determined from the partition function or from correla¬ 
tion functions. Starting from the Hamiltonian (1), the 
imaginary time path integral formalism^^ allows all ob¬ 
servables to be calculated using the partition function 
Z = J c] exp(—S'[c*, c]), with the Euclidean action 

S'[c*,c] = fg dT[c*dldT c + H]. The average number of 
particles per site, (n), is obtained by deriving with re¬ 
spect to the chemical potential, 

4 / rij exp(-S'[c*,c]) 

3 

= (Al) 

PNZ ^ ’ 
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where j is a site index and N is the total number of sites. The variation of the density with respect to the hy- 
Similarly, the compressibility is bridization yields 


d{n) _ 1 d fl d \ 

dfj, PN \Z dfj, J 


1 

' 1 d^Z 

/I azy 

PN 

Z dy? 

\z^^,) 


= dTidT2 {nj{Ti)nk{T2)) - K-(ri)) (nfe(r2)) 

jk 

= —Xq=o^uj=0, (A 2 ) 

where Xq=o,i.j„=o is to be understood as the static uni¬ 
form charge correlation function, where the limit w —> 0 
has been taken before the limit g —>■ 0. These limits do 
not commute in general^^. 


Appendix B: DMFT response formula 

In this Appendix, we show that the DMFT and DB 
response in the Hubbard model are consistent. To this 
end, we evaluate the individual terms in the expression 
(5) for the electron compressibility. 


=1/Zj D[c ,c]n—eM-S^n.,) 

- ^1/Z J D[c*,c]nexp{-Sin,p)^ 

= -l/Z J D[c*,c]ncl^c^^ exp(-S'imp) 


+ (n)'“P ( 1/Z I D[c*,c]c:^c,^exp{-S,^p) 

(B5) 


—• ^uj—0\uj—09y<y9y(y 7 


d{n) d{nr^ 

d^i dfi ^ dfi ' ^ ’ 

va 

This is done by rewriting derivatives of impurity quan¬ 
tities in terms of impurity correlation functions. Since 
the impurity problem (2) is solved (numerically) exactly, 
it is thermodynamically consistent and this rewriting is 
valid. The variation of the hybridization with respect 
to the chemical potential depends on the self-consistency 
condition (3), which includes the effect of the lattice and 
determines the derivatives. For compactness, all factors 
of P are suppressed in this derivation and the notation of 
Ref. 40 is adopted for the impurity correlation functions. 

All impurity expectation values are obtained from the 
impurity partition function 


Z = j D[c*,c]exp{-Si^p). (B2) 

For example, the impurity density is 

(n)™P = 1/Z J D[c*,c]nexp{-Simp) 

/ O 

-D[c*,c] —exp(-5'imp). (B3) 

This leads to an impurity compressibility 

= 1/Z J D[c*,c]^expi-Simp) 

= (nn)L“ro-(-)‘“"(^r'’ 

=: -Xu7=o (B4) 


where XutT,bj is the impurity three-leg vertex function^® 
and gi/cr = — the impurity Green’s func¬ 

tion. Variation of the impurity Green’s function with 
respect to the chemical potential gives 




= -llZ J D[c*,c]c^„cl^-^exp{-Simp) 

+ {^12 J D[c*,c]c^^cl^exp{-Simp)'^ 

^ y ■ 0 [c*,c]^exp(-S'ii„p)^ 

= -l/Z / I?[c*,c]c,,^c*^n(^=oexp(-5'imp) 


+ J D[c*,c\c^^cl^exp{-Simp)^ 

X ^ 1 /Z J D[c*,c]n,,;=oexp(-S'imp)^ 


— X.uj—09i'<y9i'<y 


\imp 
^1/(7,LU — O 


+ 


* ximp 

^1/(7/ 


{nY^^ 


(B6) 


The impurity Green’s function g,ya- depends on the hy¬ 
bridization even when the frequency and spin in- 














dices are different: 


dgi 


(9A, 


— =IIZ f D[c,c]c^^cl^-^ —exp(-5'imp) 

'cr' J OL\i,'(y' 

- {^12 J D[c*,c]c^^cl^exp{-Sin,p)^ 

= -l/Z J D[c*,c]c^^cl„cl,„,c^,^,exp{-Si„,p) 

+ J D[c*,c\c^^cl^exp{-Si^p)^ 

X {^IZ J D[c*,c\cl,^,c^,^,exp{-Simp)^ 

- Ir r* Ir r* \™P 

K^vg^vg) \^v' g'^v' g'I 

— '^vv' gg' ^ QvgQvgQi''g' Qy' g' 

^UI^'^GG' Ql^G^J^G-i (B7) 
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with a factor 1/N implied in the sum over k. Here 
we have introduced {^i^'a'}) for notational con¬ 

venience. The partial derivative of the self-consistency 
condition with respect to the chemical potential at fixed 
is 



= (B9) 


Here we have introduced the dual Green’s function G = 
G — g and we made use of the property 


where we have used the definition of the two-particle 
impurity vertex function (see, e.g., Ref 38). In the 
non-interacting limit, 7 = 0 and = l/{iv — 

A,,cr)- Taking the derivative —dg^J jd/S.v'a' indeed gives 

-S^„'Sacr'gi°J9^°J as in (B7). 

The DMFT self-consistency condition determines the 
hybridization function A^q. such that the local Green’s 
function is equal to the impurity Green’s function for all 
Matsubara frequencies v: 

0 — g^a ^ ^ G 

k 

^ g^a + - tk 

=U,{^X,{K'o'}) (B8) 


^ ^ ^ ^ (Givfjk T gva'){Gufj\^ -f gi/a) 

k k 

'y ^ ^i/(Tk “f 2(7i/cr y ( G^yg-k + g^a 
k k 

= G^gk + glrr- (BIO) 

k 

This holds since the DMFT self-consistency condition 
(3) implies = 0. The dependence on A^y^ re¬ 

quires some attention. The Green’s function Gjyo-k de¬ 
pends on A,yo- explicitly, but also implicitly through ^jycr. 
As shown in (B7), g^a depends on Ajy/^' also for dif¬ 
ferent Matsubara frequencies. This means that all the 
self-consistency conditions depend on all components of 
A,ycr, so the derivative of the self-consistency condition 
will be a matrix equation in frequency space: 


_df^ 

d/^y'a' 




dgva _ 

9A,,/„' ^ (n~. 


{gj + A^^-t^,y \dA^:„ 


dg, 


-1 

UG 


H" ^G 


dgi 


^ ( 1 “ 9ua X! ^ G: 


dA^.^ ^ 

-2 % 
g^a g. 


2 

i/crk 


77 ^IgIc + ^ G 


2 

i^ak* 


(Bll) 
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Collecting previous results, the variation of the self-consistency condition gives 
dfu'a' , df^'a' dA^a- 


0 = Ja1al + 

^ ^ dA,, dfx 


dA 


dfi ^ 

I/'-' 

=E 

v' g' 

=E 




-\ -1 


dA, 


dfv' 

dfj. 


_ -2 ~2 , r r , (-<2 

^ k k 


dg. 


dfj. 




“ E ^^''^'k + ^w'^aaigt'a' E *^^'cr'k 


dgi 


dg 




(mo) ^ 

v' g' 

(bj)e 


^gg' Qy' a' 


2 dgy/(ji 


dA, 


^ ^ ^l/V'k “I" SyyfScTG^gy'a' 


-■ k 


-1 


<9^1.'a- 


9^ 


“ E 


7:^'; 


7,u;—0 E ^y'cr'k + ^uv'^a 


-1 


-4 dg,icr' ^2 

ffjy'cr' / y kry/cr'l 


(B12) 


Here 


a/. 




-1 


and all other inverses should be understood as a matrix inversion in spin- and frequency space. 


Finally, the lattice compressibility is obtained by combining the above results: 
d(n) 9(n)™‘’ 9 (n)™^ i 9 Ai ,(7 


X, 


q= 0 ,cj =0 


-E 


dg. dfi ^ dAycr dg 

—Xuj—O 9" ^ ^ Xuj—0^i>G,uj—0 ^ ^ 


^2 


1 -1 


Xl^i^G^G',OJ — 0 ^ ^ ^I/ck dyyf S(jg' 

k 


(t', cj= 0 Xcj =0 


(B13) 


This is equal to the DB expression^®’^^’^®90 for y = 0 and A = 0, i.e., Xq,^ = + X^Hq^Xt,,, with dual bosonic 

self-energy 


n, 


q= 0 ,u ;=0 


k 


^i/(7k 






-1 


cr'uj—O' 


(B14) 


The ladder DB correlation function has been proven to be equal to the DMFT correlation function^®, completing the 
proof of the consistency of the DMFT response function. 


Appendix C: Second-order Dual Boson 

Eq. (B14) is the dual bosonic self-energy in the ladder 
approximation. In this paper we also use a simplier ap¬ 
proximation, which was also introduced in Ref. 40. It is 
obtained by neglecting the vertex 7 in (B14), so that 

r I 

'~fL''i‘a'{T,uj—0 ^ ^ ^u'a'k ~b ^yu'^cra' ^ ^vu' ^aa' - (Cl) 

k 

The resulting Feynman diagram (see Fig. 4) corresponds 
to a second-order approximation to the dual bosonic self¬ 
energy and reads 

TT2nd order _ \ '' \ \ '' \ (C^‘^\ 

J-iq=0,cj=0 ~ / ^ ^1^(7,cj=0 / ^ ^yGlQ^t^<7,uj—0- 

i^G k 


To get a feeling for what this approximation means, 
we can look back to the derivation of Eq. (BI4) and see 
that 7 appeared in Eq. (B7) from the dependence of g,^ 
on A,ia'- Second-order DB does not fully take this de¬ 
pendence into account, causing some degree of thermo¬ 
dynamic inconsistency. On the other hand, second-order 
DB does contain the three-leg vertices A, coming from 
both the variation of the impurity density with respect 
to the hybridization and the variation of the impurity 
Green’s function with respect to the chemical potential. 

Appendix D: EDMFT response 

As shown in (11), additional contributions to the 
electron compressibility appear in EDMET. The self- 
consistent determination of the impurity retarded inter- 
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action leads to additional terms compared to the 
result of Appendix B. Here we sketch the derivation 
and how a three- and a four-particle correlation function 
arise. The three-particle correlation function is the re¬ 
sult of the derivative /dA^^ oc — 

(n)™^ {nnY^^. The variation of the retarded interaction 
Ai^ with respect to fi is determined similarly compared 
to that of Aj^cr in (B12). The equivalent of dg^a jdA^ia' 
is, for a; 7 ^ 0, 


i.e., it gives rise to a four-particle correlator. For 
w = 0, additional terms arise, since = (nn)™^ — 
(n)™^. There are further differences in the 
derivation. The EDMFT self-consistency conditions (7) 
and (8) depend on both A^„ and via the impurity 
expectation values. This means that the formula for the 
variation of A^^ (B12) will contain contributions from 
A^„ and vice versa. Then dA^„/dg, is obtained by in¬ 
verting 


= ^ J D[c*,c]n^n^n^'n^> exp(-5'[c*, c]) 

^[c*,c]n„n^exp(-5'[c*,c]) 

X £)[c*,c]n„'ntj/ exp(-S'[c*,c]) 

= (nnnn)Z^, - {nn)Z^ (nn)LT , (Dl) 


dfv'cr' dfiyia' dAy„ dfy'a' dA^ 

dfi ^ dA^cr dg ^ dA^^ dg. 

1/(7 OJ 


(D2) 
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